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Asymptotically exact low temperature expansions for the s = 1/2 Heisenberg XXZ chains with 
boundaries are implemented by using the boundary conformal field theory. It is found that for 
1/2 < A < 1, (A, an anisotropic parameter) the boundary terms of the spin susceptibility and 
the specific heat coefficient show divergent behaviors, as temperature decreases. The degrees of 
the divergence coincide with those obtained by the Bethe ansatz method at zero temperature. 
Such low-temperature anomalous behaviors at boundaries are deeply related with finite-temperature 
corrections of the boundary entropy or the ground state degeneracy, the presence of which yields a 
very sensitive response of spin excitations. 

I. INTRODUCTION 

In the last decade, effects of open boundaries in quantum one-dimensional (ID) spin systems have been exten- 
sively studied in connection with impurity problems. 1-12 From theoretical point of view, advanced non-perturbative 
techniques such as the Bethe ansatz exact solutions, and the boundary conformal field theory have been success- 
fully applied to the investigation of boundary effects. 1-20 Although a lot of interesting features of boundary critical 
phenomena associated with open ends have been elucidated so far, there remain some issues which are not yet fully 
understood. One of them is a finite-temperature effect on the boundary spin susceptibility. According to the Bethe 
ansatz solutions for the spin- 1/2 Heisenberg chains, the SUSY t-J model, and the Hubbard model, the presence of 
boundaries gives rise to enhanced spin correlation in the vicinity of boundaries. 14 ' 15 ' 20 It was found that the boundary 
part of the uniform spin susceptibility at zero temperature for these systems behaves like ~ l/[/i{ln(/i)} 2 ] at the 
isotropic point. Here h is a magnetic field. This divergent behavior for a small magnetic field implies that near the 
boundary spin excitations are very sensitive to a perturbation induced by external fields or finite temperatures. It 
is expected that a similar divergent behavior may also appear in the temperature dependence of the boundary spin 
susceptibility. As was shown in 14,15 ' 20 , the singular contribution at zero temperature stems from the surface energy 
perturbed by leading irrelevant interactions. It was suggested in ref. 20 that at finite temperatures, in addition to the 
surface energy, a boundary entropy perturbed by irrelevant interactions also gives rise to the singular contribution. 
The boundary entropy is a residual entropy at zero temperature caused by the presence of open boundaries. 21 In the 
case that the low-energy fixed point is the Tomonaga-Luttinger liquid, namely, the Gaussian theory with the central 
charge c = 1, the boundary entropy is given by, Sb = ln(l/V2i?) for the Dirichlet boundary condition. Here R is a 
radius parameter of the Gaussian model. Thus at finite temperatures or with a magnetic field, the presence of irrel- 
evant interactions yields corrections of R and Sb , leading the singular temperature or field dependence of boundary 
quantities. 

Unfortunately, because of a technical problem inherent in boundary integrable systems, it is difficult to calculate 
boundary contributions at finite temperatures in terms of the Bethe ansatz method. However, we can compute finite- 
temperature corrections by using field theoretical methods that is based upon a low energy effective theory. For this 
purpose, the boundary conformal field theory is a very powerful tool. 22-24 In this paper, utilizing these techniques, 
we develop asymptotically exact low temperature expansions for the spin-1/2 Heisenberg XXZ chains with open 
boundaries, and calculate the leading temperature dependences of the boundary spin susceptibility and the specific 
heat coefficient. It is found that, as expected from the results obtained by the Bethe ansatz method, for 1/2 < A < 1 
(A, anisotropic parameter), these quantities show divergent behaviors, as temperature decreases. It is expected that 
the results may be relevant to experimental observations in real quasi ID spin systems with non-magnetic impurities. 

The organization of this paper is as follows. In the next section, we briefly review the low energy effective field 
theory for the spin-1/2 XXZ chain, and develop a perturbation theory in the case with open boundaries. In Sec. Ill, 
our main results for the low-temperature dependence of the spin susceptibility and the specific heat coefficient caused 
by boundary effects are presented. Discussion and summary is given in Sec. IV. 
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II. LOW ENERGY EFFECTIVE FIELD THEORY 



In this section, to establish notations, we first review the low energy effective field theory for the s = 1/2 Heisenberg 
XXZ chains briefly. 25 On the basis of this effective theory, we will perform perturbative expansions in terms of irrelevant 
interactions for the system with boundaries. The Hamiltonian for the spin-1/2 Heisenberg spin chains is given by, 

H XXZ = jJ2l S ? S i + i + S ? S ? + i + US'SUi]- (!) 

Here we consider only the antiferromagnctic case J > 0. In the massless region, < A < 1, the low energy fixed 
point of (1) is the Tomonaga-Luttingcr liquid, which belongs to the universality class of the Gaussian theory with the 
central charge c = 1. In this case, the low energy effective Hamiltonian with leading irrelevant interactions is exactly 
obtained by Lukyanov. 25 For 1/2 < A < 1, it is written as, 

H = H o + H intl (2) 

H = £ ^-[(d x <pf+U 2 }, (3) 

H mt = a 2K - 2 X [ L p cos(V8X0). (4) 
Here L is the linear system size, and the constants K, a, and A are parametrized as, 

^=[l--cos- 1 (A)]- 1 , (5) 

IT 
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The boson fields <p(x) and n(x) satisfy the canonical conjugate relation, [<j>{x), n(a/)] = iiid{x — x'). It is convenient 
to introduce the mode expansion form of 4>(x), and 0{x) = j x dx'H(x'), 
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where the operators satisfy the following commutation relations, 



(10) 



[a n ,a m ]=0, [Q,P] = [Q,P] = i. 



(11) 



N in the right-hand side of (8) is an integer corresponding to the winding number of the phase field <f>. 

As will be seen later, boundary states are constructed from the highest weight state of the U(l) Kac-Moody algebra 
\M,N), which are the eigen states of the zero mode of a n , a n defined as, 
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|M, TV) satisfies the eigen value equations, 
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a \M,N) = [<^^+X)\M,N), (13) 



2 V2KJ 

a \M,N)=(^^-^=]\M,N). (14) 



2 V2K 

The primary field corresponding to \M, N) has the conformal dimension 

r2 



ii M « + ^MN = ^ + -^r) • (15) 

In the case ofO<A<l/2 (.fT > 3/2), the low-temperature anomalous behaviors at boundaries do not appear, as 
is easily seen from the dimensional analysis. Thus we will not consider this case in the following. 

At the isotropic point A = 1 (K = 1), the effective low energy theory is described by the level k — 1 SU(2) 
Wess-Zumino-Witten model with a marginally irrelevant interaction: 

H = Hwzw + H m , (16) 

H m = -9 -^J a (*)J"(x). (17) 
Jo 2 *f-{ 

Here Hwzw is the Hamiltonian of the level k = 1 SU(2) Wess-Zumino-Witten model, and J a (x) (J a (xj) is the left 
(right) moving current of the level k = 1 SU(2) Kac-Moody algebra. The running coupling constant g depends on 
temperature T and an external magnetic field h through the scaling equation, 25 

g- 1 + \ Ha) = -Re[V(i + ^ )] + H^e^J/T), (is) 

with ip(x) the di-gamma function. 

Generally, in systems with boundaries, there may be boundary operators in addition to bulk interactions. However, 
as was pointed out in ref. 12 , in the absence of symmetry-breaking external fields at boundaries, we can exclude this 
possibility for the Heiscnberg XXZ chains. 

In the following, we will implement a perturbative expansion of the free energy in terms of leading irrelevant 
interactions, and evaluate 1 /L-corrections characterizing interesting boundary effects. For this purpose, following the 
idea of Cardy and Lewellen, 24 we consider the geometry of a semi-infinite cylinder with perimeter 1/T. Interchanging 
space and time coordinates, we define the phase field on this geometry, 

<f> c (x, t)=Q + nTPx + 2nT^L + % - V I( a „ e -^ T "(-+*) + ^ e -<w»M). (19) 

V2K 2 n 

Then, the Hamiltonian on the semi-infinite cylinder is written as, 

H C = HC + H c nt , (20) 

r 1/T dt 

= j o 2^[W) 2 + (0x* c ) 2 ], ( 21 ) 

m nt =a 2K - 2 \J^ -cos(V8K<t>). (22) 

We express the partition function by using the transfer matrix exp(—LH c ) and the boundary state \B). The lowest 
order terms of the free energy for (2) are given by, 

F _ _ £ Hole -w m + ^f dl ^(-LHS}^HS)H^ ( -,H Sm 

L L J {U\ exp{-LH )\±l) 

where |0) is the ground state of Hq. The first term of the right-hand side of (23) is the free energy of the c = 1 
Gaussian model. The second term is the 1/L correction that emerges as a result of boundary effects. Using Cardy 
and Lewellen's method, 24 we can compute this term as, 



L 



L 2na (0\B) J dX [sinh(27rTx)]^ ' (24) 

where |<&) is the primary state that corresponds to the conformal field exp(i^/8K<j)). j/$ is its conformal dimension. 
For the isotropic case, a similar consideration is applicable. In the following, we evaluate (24) exactly for particular 
boundary conditions. 
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III. FINITE-TEMPERATURE EFFECTS ON BOUNDARY CRITICAL PHENOMENA 



A. The anisotropic case 1/2 < A < 1 



In the evaluation of the boundary part of the free energy (23), we utilize properties of the boundary state. A 
conformally invariant boundary condition is imposed by demanding T(z) = T(z) at the boundary. 23 Here T(z) {T{z)) 
is the holomorophic (anti-holomorophic) part of the stress energy tensor. For the Gaussian model with c = 1, this 
condition leads the following constraint on the boundary state, 26-29 

(a„ ± a- n )\B) = 0, (25) 

where the plus (minus) sign corresponds to the Neumann (Dirichlet) boundary condition. The solution for (25) is 
expressed in terms of the Ishibashi state. 30 ' 24 For the Dirichlet condition, it is given by, 26-29 

„\ 1/4 oo , oo 



\ D )= T E e-^^exp -^i^i |M,0). (26) 

V Z ' M=-oo V n=l 71 ' 

Note that the prefactor (if/2) 1/7 4 is the boundary entropy at zero temperature. The constant 4>o is the eigen value of 
the boson field <j> at the boundary. The boundary state for the Neumann condition is, 26-29 

/ i \ !/ 4 00 / °° ~ \ 

\ N ) = (tk) E e-^^ exp ^^zIi)|o^), (27) 



2K I K^ 1 n 

' N=-oo X Tl=l 

where 9q is the boundary value of the dual field, 6(x). The prefactor (1/2K) 1 / 4 is the boundary entropy for this 
boundary condition. 

The matrix element between |0) and \B) appeared in (24) can be computed by using (26) and (27). For the case 
of 1/2 < A < 1 (1 < K < 3/2), the leading irrelevant interaction (4) is expressed in terms of the primary field 
exp(zV8if 4>) which has the conformal dimension 2K. From (15), we see that |<&) is the primary state |2, 0). (&\B) is 
non- vanishing, only if \B) contains |2, 0). The Neumann boundary state (27) does not satisfy this condition, leading 
($| N) = 0. On the other hand, ($\D) gives a finite contribution. In the following, we put (f> n — 0. This choice of the 
averaged value of the phase at the boundary is consistent with the free open boundary condition for the boundary 
magnetization. Then, we have, 

lm - (28) 

In the presence of an external magnetic field h, we add the Zeeman energy term, 

h [K rL 



h z = -s z h = --y y / dxd ^ ( 29 ) 







to the Hamiltonian (2). The free energy is evaluated by shifting the boson field 4>{x) to <f>{x) = cp(x) — yjK/2hx. 
Then the integral in (24) is replaced by, 



L 



L , cos(2Khx) 

, dx- — -i 30) 

l [smh(2TrTx)]y^ ' 

Combining (28) and (30), and carrying out the integral, we obtain corrections to the boundary part of the free energy, 

^ B = -^(2™ T ) 2K ~ 1 MB(K + 1 - 2K% (31) 

where B(x, y) = T(x)T(y)/T(x + y). 

Using (31), we can calculate the boundary contribution to the spin susceptibility, 



d 2 SF B 

Xb = 



Oh 2 
Xa 2 K 2 



h=0 

2K-3 



L 4tt 



B(K, 1 - 2^)[tt 2 - 2^'(K)](2TraT) 2K - 3 , (32) 
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where ip'(x) = dip(x)/dx. Note that for 1 < K < 3/2 (1/2 < A < 1), the boundary spin susceptibility xb shows a 
divergent behavior ~ 1/T 3 ~ 2K , as temperature decreases. This anomalous temperature dependence is also observed 
in the boundary part of the specific heat coefficient computed as, 

Cb d 2 SF B 
T dT 2 

= ^— (2K - l)(2K - 2)B{K, 1 - 2K )(2iraT) 2K - 3 . (33) 

In the above derivation, we have neglected the boundary terms of the low energy fixed point which are regular in h 
and T. We would like to stress that in the formulas (32) and (33) there is no free parameter, and the prefactors are 
exactly obtained. These divergent behaviors are physically understood as follows. In contrast to the bulk Heisenberg 
chains in which the ground state is a spin singlet state, spin singlet formation in the vicinity of boundaries is strongly 
disturbed by thermal fluctuation, because of the enhanced correlation at the boundaries which stems from the ground 
state degeneracy. It should be emphasized that the singular behaviors are not due to the presence of boundary 
operators, but interpreted as a consequence of finite-temperature corrections of the surface energy and the boundary 
entropy ln(0|_B) caused by irrelevant interactions. 

At zero temperature with a finite magnetic field, a similar singular behavior appears in the field dependence of the 
boundary spin susceptibility given by, 

X b(T = 0) = - ^ aK ] 2K ^ (2K l)r(l - 2K)h 2K -\ (34) 

L 47T 

The zero temperature susceptibility is also derived from the Bethe ansatz exact solution by using the Wiener-Hopf 
method. We have checked that (34) coincides with the result obtained by the Bethe ansatz. 



B. The isotropic case A = 1 



In the isotropic case K — 1, let us consider the boundary condition that preserves the SU(2) spin rotational 
symmetry, J a (x = 0) = J a (x = 0), which leads, 23 



{.r n + r n )\B)=Q 1 



(35) 



for a = 1, 2, 3. Here ( J°) are the n-th mode component of the left (right) moving SU(2) current. The representation 
of \B) for this boundary condition was precisely discussed by Gaberdiel and Gannon. 31 However, instead of using 
the explicit formula of \B), we exploit the SU(2) current algebra to evaluate the prefactor of (24). The marginal 
interaction (17) consists of the first descendants in the Verma module of the identity operator. The corresponding 
state vector is |$) = J1 1 J'!l 1 |0). Using the commutation relations for the level-1 SU(2) Kac-Moody algebra, 32 



[J, 



(36) 



and (0| Jt l = 0, we have, 



\JfJf\B) 
(0\B) 



\J?J-i\B) 
(0\B) 



a—l a—1 

In the evaluation of the integral in (24), we circumvent unphysical divergence by taking the limit, 



1)1^5(^,1-2^) = -1. 



Finally, we end up with, 



3aT 

U7- 



(37) 



(38) 



(39) 



In the above derivation, we take into account effects of an external magnetic field only through the field dependence 
of the coupling constant g, because it gives leading corrections induced by the magnetic field. Using (18) and (39), 
we obtain the leading term of the boundary spin susceptibility and the specific heat coefficient, 
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XB = - 



2aV/'(l) 1 
16n 2 L T[\n(aJ/T)] 



2 ' 



(40) 



C B _ 3o 1 
~T~ ~~ 4LT[ln(aJ/T)] 



2 ' 



(41) 



where ^"(1) = —2.40411..., and a = v / 27rcxp(l/4+7) with 7 the Euler constant. The temperature dependence of these 
boundary quantities is analogous to the field dependence at zero temperature, ~ l//i[ln(ft,)] 2 , which is obtained from 
the Bethe ansatz solution for the supersymmetric t-J model and the Hubbard model. '' Because of the logarithmic 
factor, the divergent behaviors of (40) and (41) arc observable only at very low temperatures T <C aJ <~ 5.73251J. 
Note that the expression (39) is applicable for h -C T. Thus we cannot derive the zero temperature susceptibility 
from (39). We need to compute higher order corrections in terms of h to take the limit T — > 0. This task is not 
straightforward because the presence of a term like ~ |ft.|/[ln(ft,)] 2 in the free energy implies that the perturbative 
expansion in h is singular. It is required to take the limit K — > 1 from the anisotropic point carefully, as was done for 
the bulk system by Lukyanov. 25 

A logarithmic singularity was also found in the NMR relaxation rate at the edge of the system calculated by Affleck 
and Qin, and Brunei et al., who showed that the bulk marginal interaction gives rise to an anomalous dimension 
at boundaries. 11,12 The above logarithmic singularity has the same origin as that for the NMR relaxation rate at 
boundaries. 



The Curie-like temperature dependence with logarithmic corrections of the boundary spin susceptibility (40) may 
be relevant to experimental observations. According to experimental measurements of the spin susceptibility for the 
Heisenberg spin chains such as Sr 2 CuC>3, a Curic-likc behavior at low temperature region is always observed, but 
has been regarded as an extrinsic impurity effect. 33 However, the results obtained in this paper imply that such an 
effect may be intrinsically important for the Heisenberg spin chains with an open end. In particular, it probes the 
character of leading irrelevant interactions in the bulk. To observe the effects obtained in this paper experimentally, 
more systematic studies distinguishing intrinsic and extrinsic effects are required. 

Although the boundary spin susceptibility at finite temperatures was studied before by both numerical and field 
theoretical methods, the anomalous temperature dependence has not been found so far. 7 ' 8 ' 10 Actually, the divergent 
behaviors of the boundary contributions (40) and (41) are visible only at very low temperatures, because of the 
logarithmic reduction factor and the small magnitude of the prefactors. It seems that these equations do not contradict 
with the numerical results obtained in refs. 8 and 10 . In the previous works based on the field theoretical methods 7 ' 8 , 
effects of irrelevant interactions on the surface energy and the boundary entropy, which are the most important 
ingredient for the divergent behaviors, have not been calculated explicitly. 

The results obtained here imply that the temperature dependences of boundary quantities similar to (40) and (41) 
should be also found in the ID Hubbard model, since its low energy effective theory for the spin sector belongs to 
the same universality class as that of the spin-1/2 isotropic Heisenberg chains. The results derived by the Bcthc 
ansatz method at zero temperature for the Hubbard model is consistent with the present results for the isotropic case, 
supporting the presence of anomalous behaviors at boundaries. 

In summary, we have obtained the boundary contributions of the spin susceptibility and the specific heat coefficient 
for the spin-1/2 Heisenberg chains using the exact effective low energy theory and the boundary conformal field 
theory. These quantities show divergent behaviors as the temperature is lowered, indicating strongly enhanced spin 
correlation in the vicinity of the open boundary. 
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